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We point out that the kinematical twist-3 contributions to the DVCS amplitude, required to 
restore electromagnetic gauge invariance of the twist-2 amplitude up to 0{t/q^), can be under- 
stood as a spin rotation applied to the twist-2 quark density matrix in the target. This allows for 
a compact representation of the twist-3 effects, as well as for a simple physical interpretation. 
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Various authors have recently studied deeply virtual 
Compton scattering (DVCS) in QCD beyond the lead- 
ing twist level . The main motivation has been the 
fact that the twist-2 contribution to the DVCS ampli- 
tude alone, originally considered in Refs. [Q-^, is not 
transverse (electromagnetically gauge invariant) even 
when neglecting terms of order t/q^ . It was shown in 
Refs. that transversality to this accuracy is re- 

stored when one includes in the light-cone expansion 
certain "kinematical" twist-3 contributions which ap- 
pear from total derivatives of twist-2 operators. Al- 
ternatively, these transversality-restoring terms can be 
represented as certain Wandzura-Wilczek type contri- 
butions to the general twist-3 skewed parton distribu- 
tions parametrizing the non-forward matrix elements 
of the relevant light-ray operators beyond the leading 
twist level 

In this Brief Report we point out that the kinematical 
twist-3 contributions to the DVCS amplitude required 
by transversality can be understood as the result of a 
certain quark spin rotation applied to the twist-2 part 
of the quark density matrix in the target. This spin ro- 
tation is necessary in order that the quark density ma- 
trix satisfy the free-field Dirac equations (up to terms 
of order t) with respect to the "incoming" and "outgo- 
ing" quark, which is both necessary and sufficient for 
the tree-level DVCS amplitude to be transverse. This 
observation allows for a very compact representation of 
the twist-3 effects within the non-local light cone ex- 
pansion, as well as for a simple physical interpretation. 
Notations and conventions in this report closely follow 
Refs. ig. 

To begin, let us recall how transversality of the tree- 
level Compton amplitude comes about in a theory of 
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FIG. 1. The virtual Compton amplitude for a free quark. 

free quarks (of charge unity, for simplicity). With the 
vector current given by J^{x) = ^p{x)^fj,ip{x), the matrix 
element of the first contraction of two currents, centered 
at point X and separated by a distance z, can be written 
in the form (see Fig.|l|) 

{p - r/2\ iTJ^iX - z/2)MX + z/2) \p + r/2)f-<= ^■ 

= itr[7^5(z)7,Af(z|X)] + {fi^u,z^-z). (1) 

Here S{z) = z/27r^(z^ — iO)^ is the free Dirac propagator 
in coordinate space (i = "/aZa), and M the density ma- 
trix constructed from the incoming and outgoing quark 
wave functions (plane waves). 



Mij{z\X) — UiUj e 
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where the spinors satisfy (p -\- f/2)u — and 
?2(p — f/2) = 0. Because of current conservation the ex- 
pression on the R.H.S. of Eq.(0) should vanish when 
contracted with the derivatives 
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FIG. 2. The twist-2 contribution to the tree-level 
hadronic DVCS amplitude in QCD. The shaded blob rep- 
resents the twist-2 part of the quark density matrix in the 
hadron, (B). 
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Since the free Dirac propagator satisfies 
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^ S{z) = -^S{z)^, = -z<5(4)(2), 



dz 
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a necessary and sufficient condition for transversality 
is that the quark density matrix satisfy the "left" and 
"right" Dirac equations 
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The free-field density matrix (g) clearly satisfies these 
equation, and thus the R.H.S. of Eq.(|l]) is transverse. 

Turning now to QCD, the twist-2 contribution to 
the tree-level light-cone expansion of the time-ordered 
product of two vector currents has a form analogous 
to Eq.(0), the only difference being that the free quark 
density matrix is replaced by the twist-2 part of the 
non- forward matrix element of the quark density tp{X — 
z/2) . . . Tp{X+z/2) between hadronic states (see Fig. ^). 
The latter is defined as (we suppress the polarization 
quantum numbers of the hadronic states for brevity) 
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where the "traceless" string operator satisfies mm 
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(7) 



and similarly for the operator with Dirac matrix £75. 
Only the chirally even part of the density matrix has 
been included in Eq.(||); the chirally odd (transversity) 
part with Dirac structure Uap would enter the DVCS 
amplitude only at power-suppressed level. Noting that 
the coefficient function for the twist-2 contribution in 
QCD is the same as that in the free theory — the free 
quark propagator, S{z) — it is evident that transversal- 
ity would again require that the quark density matrix 
satisfy the free Dirac equation. Substituting Eq. (^) in 
Eqs.(p[) one finds 
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where the ellipsis stands for the corresponding term 
with the operator ipz^^ip^ which we have not written 
out for brevity. The first term in the bracket vanishes 
because the string operator satisfies Eq.(0). The sec- 
ond term, involving the contracted derivatives, is pro- 
portional to t = or rn^ (the target mass squared) and 
can be dropped if we are interested in the DVCS am- 
plitude up to terms of order t/q^ or rv? jq^ . The third 
term, however, involving the structure 
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where ctq/j = (1/2)[7q, 7^], is not proportional to ei- 
ther t or m^, and thus cannot be neglected in DVCS 
kinematics. Thus, the twist-2 density matrix (^) does 
not satisfy the free Dirac equation even when neglect- 
ing terms proportional to t or m^. This is the reason 
why the twist-2 contribution to the DVCS amplitude 
is not transverse up to terms of order t/q^ or m? /q^ , as 
has been observed within various different approaches 
in Refs. 

One may now ask how the twist-2 part of the quark 
density matrix (^) could be modified such that it satis- 
fies the free-field Dirac equations up to terms of order 
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t or m?. Since the problematic term in Eq.(H) involves 
the Dirac matrix aap, it seems natural to look for a 
transformation in the form of a quark spin rotation. 
Since we require that the transformation matrix should 
reduce to unity for forward matrix elements, the only 
possibility (up to a variable factor in the exponent) is a 
transformation of the form 



E(z/2 
d 



exp 



1 



d 



za- 



dX 

which satisfies 

E(z/2)-i 



ZaCTajl 



E(-z/2). 



(10) 
(11) 

(12) 



It is easy to see that a density matrix with the desired 
properties is obtained by the following transformation: 
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where v = \ — v. Indeed, by a straightforward calcu- 
lation, making use of the explicit form of the matrix 
exponential 
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one can convince oneself that the "rotated" twist-2 den- 
sity matrix ([l3| ) satisfies the free-field left- and right 
Dirac equations (^ up to terms involving contracted 
derivatives of the string operator of the form 
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which have matrix elements of order t. Diagrammat- 
ically, the spin rotation of the twist-2 density matrix 



FIG. 3. Diagrammatic representation of the quark spin 
rotation, Eqs.(p^ and (^6|), as an intermediate step be- 
tween the twist-2 quark density matrix in the hadron and 
the free-quark virtual Compton amplitude. 



may be represented as in Fig.^, as a kind of "evolution 
kernel" acting between the twist-2 matrix element and 
the free quark propagator in the coefficient function. 
Note that since in the matrix element the total deriva- 
tive operator just turns into the momentum transfer, 
d/dXa = ira, the "rotated" density matrix does not 
involve any information beyond the basic twist-2 ma- 
trix elements, as parametrized by the twist-2 skewed 
distributions. 

The contribution to the light-cone expansion gener- 
ated by the "rotated" density matrix ( p^ may be writ- 
ten in the form 

hT[j^S{z)j,M{z\Xy°'] + {p^v,z^-z). 
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(16) 
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where we have substituted the expUcit form of the free 
coordinate-space quark Green function. In this expres- 
sion the spin rotation appears as acting on the coef- 
ficient functions rather than the matrix elements; this 
is of course a matter of interpretation. One can eas- 
ily verify that E g. jig ) is identical to the expression in 
Eq.(3.51) of Ref. |6|, which was obtained by making re- 
peated use of an string operator identity following from 
the QCD equations of motion and neglecting quark- 
gluon operators. Thus, we see that the "kinematical" 
twist-3 contributions to the non-forward light-cone ex- 
pansion (and consequently to the DVCS amplitude) can 
be understood as the effect of a quark spin rotation ap- 
plied to the twist-2 part of the quark density matrix. 

We have chosen here to present the spin rotation at 
the level of non-forward matrix elements of the quark 
density, as relevant to DVCS. It is clear, however, that 
this representation of "kinematical" twist-3 contribu- 
tions in the light-cone expansion is valid at the oper- 
ator level, and thus can be applied also in the context 
of other exclusive reactions where the relevant matrix 
elements are of distribution amplitude type. 
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